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The Kosterlitz-Thouless transition for the spin 1/2 Heisenberg chain with the next-to-the-nearest-neighbor 
interaction is investigated in the context of an infinite matrix product state algorithm, which is a generalization 
of the infinite time-evolving block decimation algorithm [G. Vidal, Phys. Rev. Lett. 98, 070201 (2007)] to 
accommodate both the next-to-the-nearest-neighbor interaction and spontaneous dimerization. It is found that, 
in the critical regime, the algorithm automatically leads to infinite degenerate ground-state wave functions, 
due to the finiteness of the truncation dimension. This results in pseudo symmetry spontaneous breakdown, 
as reflected in a bifurcation in the ground-state fidelity per lattice site. In addition, this allows to introduce a 
pseudo-order parameter to characterize the Kosterlitz-Thouless transition. 
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I. INTRODUCTION 

Quantum phase transitions (QPTs) 1 is one of the most in- 
triguing research subjects in condensed matter physics. A 
QPT occurs at absolute zero due to quantum fluctuations in 
a variety of quantum many-body systems. In the conven- 
tional Landau-Ginzburg-Wilson paradigm, a phase transition 
accompanies spontaneous symmetry breaking (SSB) 2 ' 3 that is 
characterized by a local order parameter. However, it now be- 
comes clear that not all QPTs fall into this category 4 . In fact, 
topological phase transitions do not involve any SSB. This 
even dates back to the Kosterlitz-Thouless (KT) transition 5 , 
first discovered for two-dimensional classical XY model. Ac- 
tually, the KT transion is ubiquitous in one-dimensional quan- 
tum systems. It describes the instability of the Luttinger liquid 
under a marginal perturbation. Normally, it is not an easy task 
to determine whether or not the KT transition occurs in a spe- 
cific system, because there are pathological problems to ana- 
lyze the KT transition numerically. One of these problems is 
that the finite size scaling technique 6 , which is successful for 
second order QPTs 1>4 , can not be applied to the KT transition 7 , 
since there are logarithmic corrections from the marginal per- 
turbation. 

Recently, a novel approach to QPTs in quantum many -body 
lattice systems has been put forward 814 , which is based on fi- 
delity, a measure of quantum state distinguishability, in quan- 
tum information science. As argued in Refsi2ii£, the ground- 
state fidelity per lattice site is able to capture quantum crit- 
icality underlying many-body physics in condensed matter. 
This fact, combining with a practical means to compute the 
ground-state fidelity per lattice site for infinite-size quantum 
lattice systems, make it practical to investigate critical phe- 
nomena in quantum many-body systems. Coincidentally, re- 
cent developments in the context of the tensor network (TN) 
algorithms for translation-invariant quantum lattice systems 
offer such a practical means. Here, we mention the infinite 
matrix product state (iMPS) algorithm 15 in one spatial dimen- 
sion and the infinite projected entangled-pair states (iPEPS) 16 
in two or higher spatial dimensions. These algorithms exploit 
the translation invariance of the system and parallelizability of 



a TN representation of quantum many-body wave functions, 
which provides an efficient way to classically simulate quan- 
tum many-body lattice systems. 

In a previous worki^, we have succeeded in applying the fi- 
delity per site approach to study the KT transition in both the 
one-dimensional spin 1/2 XXZ model and the spin 1 XXZ 
model with uniaxial single-ion anisotropy. This results in 
the introduction of a novel concept-/weMiio SSB, offering a 
novel perspective to understand the KT transition, in the con- 
ventional Landau-Ginzburg-Wilson paradigm, from the iMPS 
representation. As such, it resolved the controversy regarding 
whether or not the ground-state fidelity is able to detect the KT 
QPTs 13 . However, more work is needed to clarify if such an 
approach is applicable to the KT transition in other quantum 
many-body lattice systems. 




FIG. 1 : (Color online) A sketch for the zigzag chain with the nearest- 
neighbor (NN) and next-to-the nearest-neighbor (NNN) interactions 
J and J<5. The parameter 6 represents the ratio between the NNN 
coupling and the NN coupling. The dash-line box indicates the unit 
cell. 

In this paper, we propose an iMPS algorithm, a general- 
ization of the infinite time-evolving block decimation algo- 
rithm 15 , which allows us to accommodate both the next-to- 
the-nearest-neighbor interaction and spontaneous dimeriza- 
tion. It is found that, in the critical regime, the algorithm 
automatically leads to infinite degenerate ground-state wave 
functions, due to the finiteness of the truncation dimension. 
This results in pseudo symmetry spontaneous breakdown, as 
reflected in a bifurcation in the ground-state fidelity per lat- 
tice site. In addition, this allows to introduce a pseudo-order 
parameter to characterize the KT transition, which must be 
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scaled down to zero in order to be consistent with the Mermin- 
Wagner theorem^. 



II. MATRIX PRODUCT STATE ALGORITHM ON AN 
INFINITE-SIZE ONE-DIMENSIONAL LATTICE 

Suppose the model Hamiltonian takes the form: H = 
Yjih [u ' +2] , with h v ' i+2] being the sum of the nearest- 
neighbor and the next-to-the-nearest-neighbor three-body 
Hamiltonian density. A realization of such a model Hamil- 
tonian in the zigzag chain is shown in Fig. Q] To take ac- 
count of a possible dimerization arising from the competition 
between quantum fluctuations and geometric frustration, we 
choose four sequential sites (A, B, C, D) as a unit cell. Start- 
ing with a randomly chosen initial state ^(O)), which is not 
orthogonal to the genuine ground state, a ground-state wave 
function can be computed by the imaginary time evolution 
Wg) = exp(-//r)| , F(0)>/||exp(-//r)| , F(0))|| with r -> oo. To 
realize the imaginary time evolution operation, the imaginary 
time t is divided into many small slices 6t = t/N to ap- 
proximate the continuous time evolution by a sequence of 
small gates. Meanwhile, the time evolution operator is ex- 
panded to a product of evolution operators acting on three 
sites: Uj = exp(-/i [ '' ,+21 5r) with 6t « 1, as follows from 
the Suzuki-Trotter decomposition 19 . Any wave function ad- 
mits an iMPS representation in a canonical form: attached to 
each site a four-index tensor Y^ r and each bond a diagonal 
matrix A„, where n — A,B,C and D. Here, s is a physical 
index, s = 1, • • • , d, with d being the dimension of the local 
Hilbert space. / and r denote the bond indices, /, r — 1, • • • ,x, 
with x being the truncation dimension. For simplicity, we use 
r„ instead of r_ in the text bellow. 
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FIG. 2: (Color online) The updating procedure for the tensors T„ and 
the diagonal matrices A„ (n = A,B,C and D): (i) Apply the three- 
site gate f/3. one needs to update the tensors T A , T B , Tc, and the 
diagonal matrices A A , A B ; (ii) Contract the tensors inside the dash- 
line box in (i) into a a single tensor M; (iii) Reshape M into a matrix 
and perform a singular value decomposition (SVD); (iv) Insert the 
identity resolution on the right hand side; (v) Update the tensor Tc', 
(vi) Perform an SVD for the matrix contracted from M\ and A B \ (vii) 
Insert the identity resolution; (viii) All the tensors T A , T B , Tc, and the 
diagonal matrices A A and A B are updated. 



The updating procedure for the tensors F„ and the diago- 
nal matrices A„ (n = A,B,C and D) in the iMPS representa- 
tion under the action of the three-site gate f/3 is visualized in 
Fig- 12 (i) Apply the three-site gate f/3 onto the iMPS tensors 
Fa, Tb and Tc- The tensors involved in this action is shown 
in the dash-line box. (ii) Contract the tensors Aq, F^, Aa, r#, 
Ab, Fc and Ac into a single tensor and reshape the tensor into 
a^d 2 X%d matrix M. (iii) Perform a singular value decompo- 
sition (SVD) to the matrix M. After truncating and reshaping, 
we get the tensors M\, V\ and the updated diagonal matrix 
Ab- (iv) Insert the identity resolution / = A c l Ac on the right 
hand side and contract the tensors V\ and A c l to get the up- 
dated fc, as done in step (v). (vi) Contract M\ and Ab into a 
matrix and perform a SVD again, one gets the tensors V2, V3 
and the updated Aa- (vii) Insert the identity resolution; (viii) 
All the tensors P4, f b, fc, and the diagonal matrices Aa and 
Ab are updated. Shifting the action of the three-site gate t/3 
one site each time and repeating four times, we are able to up- 
date the tensors under the imaginary time evolution of an slice 
St. Repeat the procedure until the ground-state energy per lat- 
tice converges, an approximate ground-state wave function is 
generated in the iMPS representation. 



III. MODEL 

Consider a spin 1/2 Heisenberg chain with the nearest- 
neighbor coupling J and next-to-the-nearest-neighbor cou- 
pling J6. It is described by the Hamiltonian 
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where S w are the spin-1/2 Pauli operators at the 2-th site. The 
system is equivalent to a zigzag chain as shown in Fig.Q] We 
set the nearest-neighbor antiferromagnetic coupling J — 1 as 
the energy scale and consider the ratio interval < 6 < 0.5. 
The system undergoes the KT transition at 6 C ~ 0.241 1—: the 
Luttinger liquid state for 5 < 6 C and the dimerized state for 6 > 
6^, respectively. Actually, the KT transition accompanies a 
discrete Z2 SSB in the dimerized phase. 



IV. GROUND-STATE FIDELITY PER LATTICE SITE 

The ground-state fidelity per lattice site, d(5\ , 62), with 6 be- 
ing the control parameter, is defined as the scaling parameter: 
F{6u&x) ~ d(6i,6 2 ) N , with N the total number of the lattice 
sites, and F{6\,62) = {(^(Si^iSi)^ between two ground- 
state wave functions ^(^i)) and ^(62))- It characterizes how 
fast the fidelity F (6\ , 62) varies when the thermodynamic limit 
is approached 9 ^. In fact, the ground-state fidelity per lattice 
site d(S\ , 62) may be regarded as a partition function per site of 
a classical statistical vertex lattice model defined on the same 
lattice. That explains why d(6\,d2) is able to detect QPTs 22 . 
That is, the ground-state fidelity per lattice site exhibits a sin- 
gular behavior when the control parameter 6 crosses a transi- 
tion point 6 C . 




FIG. 3: (Color online) The ground-state fidelity per lattice site, 
d(6i,52), for the spin 1/2 Heisenberg chain with the NNN interac- 
tion. We have chosen ^(62) as a reference state, with 62 in different 
phases: (a) <5 2 = in the S 1/(2) symmetry-broken phase. There 
is a bifurcation point in d{8\,8>i), which tends to disappear, when^ 
approaches 00. (b) 62 = 0.4 in the Z2 symmetric-broken phase. A 
bifurcation point always exists whatever the truncation dimension x 
we choose, as it should be for a discrete symmetry SSB. Therefore, 
the ground-state fidelity per lattice site, d((5i,<5 2 ), is able to distin- 
guish degenerate ground states, with a (pseudo) critical point as a 
bifurcation point. 



Fig. [3] shows the ground-state fidelity per lattice site, 
d(5i, 62), as a function of Si, when 62 is fixed, for the Heisen- 
berg model with the NNN interaction, for different values 
of the truncation dimension x- The iMPS simulation is per- 
formed for a randomly chosen initial state. It automatically 
induces degenerate ground-state wave functions, which break 
the S U(2) symmetry in the Luttinger liquid phase and the Z2 
symmetry in the dimer phase. 

In Fig. Oa), the reference state X ¥(S2 = 0) is chosen in 
the S U(2) symmetry -broken phase. The ground-state fidelity 
per lattice site, flf(c>i, £2), takes infinitely many values lying 
between two extreme values, as a consequence of degener- 
ate ground states, whereas in the Z2 symmetry-broken phase, 
d{6\,62) yields just one value. A bifurcation point occurs, 
which coincides with the pseudo-transition point With 
increasing x, the difference between two extreme values of 
d{5\,62) decreases. Note that, the bifurcation point tends 
to disappear when x g° es to infinity, implying that degener- 



FIG. 4: (Color online) The pseudo-local-order parameter O x as a 
function of <5. The iMPS simulation is performed for a randomly 
chosen initial state. An extrapolation of the pseudo critical point S x c 
is performed for pseudo-transition points, at which the pseudo-order 
parameter become zero, yielding the KT transition point S c = 0.2418 
in the right-top inset. The right-down inset shows that, in the critical 
phase, the pseudo-order parameter is scaled down to zero according 
to O x = ax~ b 0- + c x ~ 1 \ with a=0.5993, b=0.6150, c=1.0379, to keep 
consistent with the Mermin-Wagner theorem. 



ate ground states arise from the finiteness of the truncation 
dimension, an artifact of the iMPS algorithm. In fact, all 
the degenerate ground states should collapse into the genuine 
ground state as x approaches 00, as required to keep consis- 
tent with the Mermin-Wagner theorem: no continuous sym- 
metry is spontaneously broken in one-dimensional quantum 
systems. However, the critical point 6 C may be determined by 
performing an extrapolation with respect to a few reasonably 
small x's. 

In Fig. [3jb), the reference state ¥(62 = 0.4) is chosen in the 
Z2 symmetry -broken phase. Bifurcation points occur again in 
the ground-state fidelity per lattice site. However, they never 
vanish, when x approaches 00. Instead, they tend to saturate 
when x increases. This is expected, since bifurcation points 
arise from the discrete group Z2 SSB. 

Therefore, we conclude that, for a finite truncation dimen- 
sion x> the ground-state fidelity per lattice site, d{5\,62), is 
able to distinguish degenerate ground-state wave functions 
with a reference state in the symmetry-broken phase. In ad- 
dition, the iMPS algorithm enables us to locate the KT tran- 
sition point accurately by computing the ground-state fidelity 
per lattice site, d{5\,62), with a moderate computational cost. 



V. PSEUDO ORDER PARAMETER 

As is well known, the iMPS algorithm yields the best ap- 
proximation to a ground-state wave function for a gapful lat- 
tice system. However, for a gapless system with continuous 
symmetry, if the truncation dimension x is finite, the iMPS 
algorithm automatically produces infinitely many degenerate 
ground states from a randomly chosen initial state 17 , each of 
which breaks the continuous symmetry. That is, a numeri- 
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cal phenomenon occurs, which shares all the features of an 
SSB resulting from random perturbations. Such a pseudo- 
symmetry-broken order may be quantified by introducing a 
pseudo-local-order parameter, which may be read off from a 
reduced density matrix on a local area^. However, this phe- 
nomenon is in contradiction with the Mermin-Wagner the- 
orem, which states that no continuous symmetry is sponta- 
neously broken for quantum systems in one spatial dimen- 
sion^. To resolve this apparent contradiction, one has to re- 
quire that the pseudo-local-order parameter must be scaled 
down to zero, when the truncation dimension ;f goes to oo. We 
emphasize that both pseudo SSB and pseMt/o-local-order pa- 
rameter arise from the artifact of the iMPS algorithm, in sharp 
contrast with a genuine SSB and a local order parameter. 

In the Heisenberg model with the NNN interaction, 
the pseudo-local-order parameter may be chosen as O x — 

■yJiS'j-) 2 + (S' y ) 2 + (S' z ) 2 for the truncation dimension x- The 
pseudo-local-order parameter O x is plotted as a function of 
6 in Fig. [4] The iMPS simulation is performed for a ran- 
domly chosen initial state, with the truncation dimension x 
to be 4, 8, 16, 32, and 64, respectively. Notice that O x is zero 
in the Z2 symmetry-broken phase, but nonzero in the pseudo 
S U (2) symmetry-broken phase, due to the finiteness of the 
truncation dimension X- However, this is nothing but an ar- 
tifact of the iMPS algorithm. Remarkably, one may take ad- 
vantage of this artifact to locate a critical point Note that 
an extrapolation with respect to x may be performed for the 
pseudo-phase-transition points df,, at which the pseudo-local- 
order parameter becomes zero. In the right-top inset of Fig. |4] 
the KT transition point 6 C — 0.2418 is determined from such 
an extrapolation, which is comparable with 6 C = 0.241 1 from 



the level spectroscopy analysis — . Here, we require that the 
pseudo-local-order parameter must vanish, when x goes to 00, 
to keep consistent with the Mermin-Wagner theorem. There- 
fore, a fitting function O x = ax (J + cx~ l ) is chosen, with 
a=0.5993, b=0.6150, c=1.0379, as shown in the right-down 
inset in Fig.@] This ensures that, when^- — » 00, all the degen- 
erate ground states collapse into the genuine ground state. 



VI. CONCLUSION 

We have investigated the KT transition for the spin 1/2 
Heisenberg chain with the NNN interaction in the context 
of an iMPS algorithm, a generalization of the infinite time- 
evolving block decimation algorithm to accommodate both 
the next-nearest-neighbor interaction and spontaneous dimer- 
ization. It is demonstrated that, in the critical regime, the 
algorithm automatically leads to infinitely many degenerate 
ground-state wave functions, due to the finiteness of the trun- 
cation dimension. This results in pseudo symmetry sponta- 
neous breakdown, as reflected in a bifurcation in the ground- 
state fidelity per lattice site. In addition, this allows to intro- 
duce a pseudo-local-order parameter to characterize the KT 
transition, which scales down to zero when x approaches in- 
finity. 
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